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Constraining the dark energy equation of state with 
double source plane strong lenses 
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ABSTRACT 

We investigate the possibility of constraining the dark energy equation of state 
by measuring the ratio of Einstein radii in a strong gravitational lens system with 
two source planes. This quantity is independent of the Hubble parameter and directly 
measures the growth of angular diameter distances as a function of redshift. We in- 
vestigate the prospects for a single double source plane system and for a forecast 
population of systems discovered by re-observing a population of single source lenses 
already known from a photometrically selected catalogue such as CASSOWARY or 
from a spectroscopically selected catalogue such as SLACS. We find that constraints 
comparable to current data-sets {cr{w) ~ 15%) are possible with a handful of double 
source plane systems. We also find that the method's degeneracy between H,m and w 
is almost orthogonal to that of CMB and BAO measurements, making this method 
highly complimentary to current probes. 
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1 INTRODUCTION 

We appear to live in a n expanding universe, and the ex- 
pansi on is accelerating (|Riess et al.l 1 19981 ; IPerlmutter et al.l 
ll999T ). This acceleration cannot be caused by the baryonic 
matter that we interact with everyday, because it exerts a 
braking effect on the universe's expansion; instead, the ac- 
celeration must be caused by something (frequently termed 
dark energy) that exerts a negative pressure on the universe. 
Our current understanding of dark energy is very poor; from 
observations we know that it to day makes up arou nd 70% of 
our universe's energy density IjLarson et a.1.1 120111 ) and has 
an equation of state with to « -1, where w is defined as 
the ratio of pressure to energy density, w — p/pc 2 . Any 
equation of state with w ^ — ^ will exert sufficient nega- 
tive pressure to cause accelerated expansion, so why do we 
currently observe w ~ — 1? 

Little is known about dark energy so there are 
many possible explanations that we cannot currently ex- 
clude. The standard model of a cosmological constant 
(with w = — 1) is consisten t with the current 

best constraints (lEisenstein et al.ll2 005: Pcrc ival et al.ll2010l ; 
iKomatsu et~aT1 1201 ll ; ISuvu et all l2010h . but thecurrent 
uncertainties leave much room for manoeuvre, includ- 
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ing models where the equation of s tate varies with time 
IjCaldwell. Dave, fc Steinhardtl ll998T ). Tighter constraints 
will allow us to better pin down the nature of dark energy 
and the role it plays in the expansion of our universe. 



First discovered 

IWalsh. Carswell. fc Wevmannl 
itational lens ing is now 



(c 



Q0957+561 by 

1 19791 ) , strong grav- 
pow erful co smological 



Witt. Mao, fc Keetonl |2000|; iKochanok 



Saha & Williams 200 j; ISchechterl 120051; lOguri 



20ld ; 



tool 
2002; 

| 2007l; ICongdon. Keeton fc Nordgrenl 12008 
Keeton fc Moustakasl 120091 ; ISuvu et all |2010T ) and the 
time delay strong lensi ng con straints on dark energy are 
forecast to improve ((Coe_& Moustakas 2009). In this 
this work we examine how well w can be constrained 
using a population of gravitational lens systems with 
two background sources at different redshifts (schematic 
shown in Figure [TJ. Such systems are known with cluster 
scale lenses fe.g.. ISoucail. Kneib. fc Golsd 12004 ; I Julio et al.1 
l2010h and a g alaxy-scale double s ource plane system was 
recently found (JGavazzi et al.ll2008l ). although cosmography 
is currently imprecise using the galaxy-scale system, as no 
spectros copic redshift for the fa r source has been obtained 
(but see ISonnenfeld et al.1 l201ll ) . The incentive to use such 
systems is strong; current methods of co nstraining w are 
degenerate with the Hubble constant (lEisenstei n et al 



2005; Pcr cival et aLll2010l ; IKomatsu et al.ll201ll ; ISuvu et al 
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Figure 1. Sketch of a double source plane lens system. For a 
singular isothermal sphere, the dimensionless number r\ is the 
product of -Di s i and D a 2 (both in red) divided by the product of 
As2 and D B i (both in blue). 



l20ld ; IWeinberg et all |2012| ), but taking the ratio of the 
two Einstein radii in a double source system allows us 
to constrain w independently of the Hubble constant. 
This independent constraint will not only teach us about 
the nature of dark energy, but will also help improve 
constraints on other cosmological parameters by providing 
complimentary probes with a prior on w. For example the 
WMAP7 constraints on h (the reduced Hubble constant) 
are h = 0.710 ± 0.025 (|Komatsu et all l201ll ) for a flat 
ACDM model with w fixed at —1, whilst allowing w to be 
a free parameter loosens the WMAP-only constraints by 
nearly an order of magnitude to h = 0.751q'i4- 

This paper is meant to investigate the cosmological in- 
formation content of double source plane lenses; we intend 
to address the following questions: 

(i) Can double source plane lenses be used for cosmogra- 
phy? 

(ii) What is an optimal configuration of lens and source 
redshifts for cosmography? 

(iii) How well could cosmography be constrained with a 
realistic population of double source plane lenses? 



In future papers we will investigate the observational diffi- 
culties and potential systematics, as well as producing de- 
tailed forecasts for a double source plane search with up- 
coming facilities. 

The paper is laid out as follows. In Section 2 we outline 
the method and relevant quantities. Section 3 looks at the 
prospects for cosmography with a single hypothetical lensed 
double source system and should provide the reader with 
an intuition of what is an optimal system for cosmography, 
whilst Section 4 investigates cosmography with a popula- 
tion of double source plane lenses. The population of dou- 
ble sources are forecast for re-observing a pre-existing lens- 
source catalogue. In Section[5]we examine the cosmographic 
constraints for an evolving equation of state. Section 6 con- 
cludes by examining the effectiveness of combining a set of 
double source plane systems with current constraints upon 
w from other probes. 

Throughout this work we assume a spatially flat fiducial 
cosmology with w — — 1 at all times and £7m = 0.27; where 
necessary we take h — 0.7. All models assert flatness and 
have a constant equation of state, except in Section[5] where 
we investigate an evolving model and spatial curvature. 



2 GRAVITATIONAL LENSING AND DOUBLE 
SOURCE SYSTEMS 

Mass causes the deflection of light and so can act as a gravi- 
tational lens, producing a set of images where the time delay 
function is stationary. For a strong lens system with a source 
lying directly on the optical axis this produces an Einstein 
ring. The Einstein radius is sensitive to the angular diameter 
distances between observer, lens and source. In the case of a 
thin lens in an otherwise homogeneous universe the Einstein 
radius is given by 



Be = 



4GM{0 E ) As 
c 2 D ol D 



(1) 



where 9e is the Einstein radius, M(6e) is the projected mass 
within the Einstein radius, and the Du are the angular di- 
ameter distances between observer (o)ljlens (1) and source 
(s). In this work we assume all lenses to be Singular Isother- 
mal Spheres (SISs), with a spherically symmetric mass dis- 
tribution. The density of an SIS profile is given by: 



p(r) 



<Ty_ 

2nGr 2 



(2) 



SIS mass distributions have been shown to approximate 
the observational data well (e.g. iKoopmans et al.l 120091 ; 
lAuger et al.l [2010J) and allow for simple lensing forecasts 
since the lensing deflection angle and magnification have an- 
alytical forms. The Einstein radius of an SIS is given by 



sis - °V As 
9e =47r ~A 



(T V 



2 jj ls 
I -=— arcseconds. (3) 

i-'s 



V186 kms 

For an SIS model, strong lensing occurs when the source is 
in the region < 8 B < 9e, where 9 B is the unlensed angular 
position of the source with respect to the optical axis. Two 
images are produced at #+,- = 8 S ± #e with magnification 
^+,_ = l±0 B /0,. 

Since the equation of state governs the expansion of 
our universe it affects the evolution of angular diameter dis- 
tances as a function of redshift. 



D- - C/H ° 
13 (1 + 3) 




(4) 



where sinn(x) = sin(x), x, or sinh(a;) for open, flat, or 
closed universes respectively, and E(z) is the normalised 
Hubble parameter: 



E(z) = 



H{z) 
H a 



(5a) 



E wCUM = ^n M (l + zy i + n k (l + z)i + (n dc )(l + z)^+™) 

(5b) 

E W * CDM = ^/n M (l + z)3 + n k (l + zY + (fi dc )e'(*>. (5c) 

I(z) is an integral given by 

(1 + *)- (6) 

Equation I5bl holds if w is constant, whilst Equation I5cl is 
general for any universe with a time-evolving equation of 

1 From now on we drop the 'o' referring to the observer in 
observer-lens/source distances 



/(*) S3 / 1 -^^Az' 
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Figure 2. The w and J!m plane: black are the con- 
straints possible from our fiducial single double source system 
[z\ = 0.35, z s i = 0.6, z s 2 = f.5), given uniform priors 
of ^ Qm ^ 1 an d ~ 2 ^ w ^ — A and a measurement with 
cr{rj)/{rf) = 1%. Grey shows the WMAP7 results. Contours are 
1, 2 and 3 <r. The orthogonality (and hence complementarity) be- 
tween WMAP and the double source plane constraint is striking. 



state. Neglecting the mass of the closer source, we define 
the quantity rj as the ratio of the two Einstein radii, 



1 = 



(7) 



with si and s2 referring to the near and far source, respec- 
tively. For an SIS lens rj is given by 

^S = AsrA. 

Ms2i^s 1 

The ratio, rj, has the intrinsic advantage that it is indepen- 
dent of the Hubble constant and is only weakly dependent 
upon the mass distribution of the lens (in the case of an SIS 
model r\ is independent of the mass), so is a function only 
of w, Qm an d the redshifts of the lens and sources (and the 
mass model). 



3 SINGLE SYSTEM CONSTRAINTS ON A 
CONSTANT EQUATION OF STATE 

In order to investigate [w, Qm] we first fix our lens 
and source redshifts to gain intuition about the problem. 
We choose a hypothetical double source plane lens, with 
\z\ — 0.35, z s i = 0.6, z B 2 — 1.5] as the "fiducial system". 
We use a uniform prior of < Qm ^ 1 and — 2 < w < — | 
and assume the observational precision on r\ to be 1% (i.e. 
a ( r l)/{ r l) = 0.01) to perform a Markov Chain Monte Carlo 
(MCMC) analysis of the w-il M plane. For our MCMC anal- 
ysis we assume Gaussian likelihoods given by 

1 (J](Qm, W) — ?7Fiducial) 2 

'2 <J 2 
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Figure 3. The forecast 1, 2 and 3 cr confidence re- 
gion on w possible from a single double source system 
(z\ = 0.35, z s i = 0.6, z s 2 = 1.5), as a function of the un- 
certainty on r], marginalised over all fJu using a WMAP prior. 
The dotted vertical lines represent the la constraint on w from 
WMAP alone. Uncertainties significantly below 1% are currently 
unrealistic given the need to correctly model the compound lens- 
ing i ntroduced by t he fir st source. The inferred uncertainty on 
r) in iGavazzi et al.l l|200a ) is 1.2%, without accounting for lens- 
ing by the intermediate source but only using a conjugate points 
analysis rather than all available pixels from the HST images. 



where ?7Fiduciai is the value of rj expected for the fiducial cos- 
mology. The resulting 1, 2, and 3 a contours are shown in 
Figure [2] Given the assumed 1% measurement, we find that 
it would be possible to constrain w < —0.54 at the 95% CL, 
for the fiducial system. However the posterior gives little 
constraint on w < —0.7 without a tighter prior on Um; this 
occurs because in a high £7m universe, dark energy plays very 
little role in the expansion history and so rj becomes broadly 
insensitive to the equation of state. A significa nt body of ev- 
idenc e now exists to suggest Qm ~ 0.3 (e.g., lLarson et al.l 
120111 ) ; this allows us to exclude the right of Figure [2] and 
places our universe in the regime where rj is sensitive to 
the dark energy equation of state. In order to take advan- 
tage of th is prior knowledge we use the WMAP 7 year re- 
sults from lKomatsu et al.l i|201ll ) and resample their MCMC 
ch ainjj using the i mport ance sampling method as outlined 
in iLewis fc Bridle! (|2002h . We use the MCMC chains pro- 
duced assuming a flat wCDM model, giving a set of points 
each possessing a value for Hm, W and a weight. Each point 
in the WMAP chain is then given a likelihood Li, accord- 
ing to Equation [3] and the chain is re- weighted such that 
the posterior is the product of the original weight and the 
likelihood. In Figure [3] we plot the 1, 2, and 3 a confidence 
intervals for the fiducial system plus WMAP as a function of 
percentage uncertainty on 77; an uncertainty of less than 2% 
is key to improving the constraint on w with a single system. 
Whilst minimising the uncertainty on r\ maximises the con- 
straint on w, a measurement much more accurate than 1% 
is unlikely with current resources, given the need for space 
telescope time in order to constrain well the lens mass model 
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— [W-O.35, 3,1=0.6, *rf=1.5] 

— [ai-0.35, z,i=0.6, z,s=2.0] 



— [zi-0.35, 1,1=0.6, r,2=1.5] 
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Figure 4. Same as Figure [2] but with a different double source plane system over- plotted in red. Left: The second source moved further 
away (z\ = 0.35, z a i = 0.6, z s 2 = 2.0). Centre: The first source moved closer to the lens (z\ = 0.35, z s \ = 0.45, z s 2 = 1.5). 
Right: The lens and first source moved to a redshift of 0.2 greater (z\ = 0.55, z s i = 0.8, z s 2 = 1.5). One can infer from these that 
the optimal system has a close pair of lens and first source, with a distant background source. Whilst moving the lens and first source to 
higher redshift alters the tilt of the banana, there is only marginal improvement on the constraints in the region that is also consistent 
with WMAP7. 
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Figure 5. Contour plots of the symmeterised 1<t constraints on w possible for a single double source plane system combined with WMAP, 
for fixed lens position and variable source positions, given an uncertainty of a-(rj)/r] = 1%. The dashed line indicates z s i = z s 2- Left: 
Lens redshift Z\ = 0.35. Right Lens redshift Z\ = 0.55. The contours show that maximising the separation between first and second source 
and minimising the separation between lens and first source gives the best constraints on w. 



and the image separ ation. The inferred uncertainty on r\ in 
iGavazzi et al.l (|2008l ) is 1.2%, neglecting any uncertainties 
induced by lensin g from the firs t sourc e. These systems are 
compound lenses; IGavazzi et al.l i|2008l) use the Tully- Fisher 
relation to estimate the mass of the first source and con- 
servatively estimate that the uncertainty on r\ grows to ap- 
proximately 6 per cent. Their estimate does not attempt to 
model the relensing and hence neglects most of the informa- 
tion available in the high-resolution imaging, and modeling 
the lensing by the first source of a typical compound lens 
is likely to significan tly improve upon the 6 % estimate of 
Gavaz zi et al.l (J2008). It has been shown (e.g. IVegetti et al.l 
2010 fc 20121 ) that using a pixelated image modelling tech- 
nique can quantify perturbations due to substructure with 
masses as low as 0.1% of the the primary lens mass (enclosed 
within the Einstein radius), although this method has not 
yet been applied to significant perturbers outside the pri- 
mary lens plane. In this work we take cr{rj)/{rf) = 1% un- 
less otherwise stated; this may prove optimistic, but serves 



to illustrate the cosmological information encoded in such 
systems. We do not wish to imply that 1% is definitely 
achievable and we will explore the precision with which r\ 
can be measured in future work (Collett et al. In prep). 
A measurement of the fiducial system at cr{rf)/{rj) = 1% 
gives the constraint of W — — 0.991q22 when combined 
with the WMAP7 dataset. The main improvement upon the 
WMAP data comes from the fact that low w, low Qm uni- 
verses (bottom left of Figure [5J) are strongly excluded by the 
double source plane system given our fiducial cosmology. 

In order to further build intuition and to suggest an 
observing strategy, we investigate the optimal [z\, z s i, z S 2] 
configuration for a single system. Figures [4] and [5] show the 
effect of altering the redshifts of each object in the system. 
We find that increasing the redshift of the second source im- 
proves the constraint on w; this makes physical sense since 
it means D S 2 and D\ S 2 probe a longer part of the cosmic 
expansion history, amplifying any deviations from D™^~ . 
We also find that the constraint on w is improved by de- 



© 0000 RAS, MNRAS 000, [Q-?? 



Constrainina wnw. with Comvound lenses 5 



creasing the separation of lens and first source. Small Azi_ s i 
and large Az s i. S 2 being the optimal configuration is physi- 
cally sensible - it makes Di&J D\ B i most sensitive to w, since 
very little of the D\ B 2 integral is spent in the redshift range 
also included in D\ B \. It is also important that Z\ > 0.2 or 
the ratio of D B /D\ B quickly saturates to I independent of 
the underlying cosmography, ensuring r\ ~ I; this is also the 
case when A« 8 i- S 2 is small. 

Small z s i — z\ also increases the importance of lensing 
by the intermediate object and could introduce a significant 
systematic error if the lensing by the first source is modeled 
incorrectly. We leave the analysis of secondary lensing for a 
later work, but note that it should be possible to model for 
bright systems due to the wealth of positional information 
contained in two bright arcs. For low mass first sources sec- 
ondary lensing is likely to be dominated by external conver- 
gence be_tween_the_sources, which can be statistically mod- 
elled (jSuvu et al.ll2010P ). Secondary lensing is likely to cause 
problems in systems where the intermediate source is of in- 
termediate mass and the background source is faint. A sec- 
ond factor that may trouble systems with small z B \ — z\ is 
that the mass density profile within the inner Einstein ring 
may be significantly different to that contained by the outer 
ring and this will mean the uncertainty on r\ could be signif- 
icantly affected by the uncertainty on the mass distribution 
of the lens. Hence, for very small Azi_ s i it may be impos- 
sible to achieve a 1% uncertainty on 77, although the cross 
section for very small Azi_ s i is also very small. 



4 PROSPECT FOR COSMOGRAPHY WITH N 
DOUBLE SOURCE SYSTEMS 

Angular diameter distance is an integral of w and Qm over 
redshift, so the exact nature of the degeneracy between the 
two parameters varies with z\, z s i, and z B 2- The curve of the 
"banana" shaped distribution in Figure [2] depends on the 
redshifts, as does the constraining power of each system; 
hence we do not necessarily expect n systems to provide a 
i/n improvement. Indeed, the precision could scale better or 
worse than y/n depending on the widths and orientations of 
the constraints from each system; the constraining power of 
n systems is not trivial to predict. 

We now draw a set of double source systems from a 
realistic population of lenses and sources. We assume that 
double lenses are more likely to be found by observing 
already-known single lens systems to greater depth and at 
different wavelengths. We focus on two types of initial sur- 
veys: photometrically selected wide arc catalogues where the 
survey looks for multi ple sources around a central object 
(e.g., CASSOWARY: feelokurov et alJl2007h and spectro- 
scopically selected catalogues where the survey looks at the 
spectra of galaxies for higher redshift emission lines, such as 
SLACS IJBolton et al.ll2004h . 

Follow up of an existing catalogue introduces two selec- 
tion functions, the selection function of the original survey 
and of the follow up. We estimate the selection function for 
current surveys in two ways: we create a mock photometric 
catalogue for CASSOWARY-like surveys and we draw ran- 
dom s ystems from the current SLACS sample of lAueer et al.l 
(2009). For the follow up selection function we assume the 
observing strategy will prioritise those systems which are 
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Figure 6. Example probability distribution function for finding a 
second source behind a lens at z\ = 0.2, with velocity dispersion 
350 km s —1 . given observations at 1.1 mm down to a flux limit 
0.3 mjy. The PDF does not integrate to one, since it is not certain 
that a second source is found. The dashed line shows the unlensed 
number density of sources with S > 0.3 mjy per squar e degree 
(scale on right hand y-axis) in the Be thermin et al.l 1 1201 II ) model. 
Some of the S < 0.3 mjy objects also contribute to the lensed 
PDF due to magnification. 



best for cosmography. In practice this usually means look- 
ing for high-redshift background second sources given a low 
redshift (z\, Zsi < 1.5) lens-source pair, or low redshift inter- 
mediate sources, given a well separated lens and source. As 
such we assume follow up will either be conducted in the mil- 
limetre, in order to preferentially select higher redshift sec- 
ond sources, or in deep optical imaging, where low-redshift 
intermediate sources could be found. For the millimetr e pop- 
ulatio n we adopt the unlensed model of Bethermin et al. 
(|2011T ) . assuming observations at 1.1 mm and continuum 
flux limits of 0.3 mjy and 1 mjy; an example source red- 
shift probability distribution function for a lens at Z\ — 0.4 
is shown in Figure [6] For the optical source distribution 
we use the photometric redshift catalogue of galaxies ob- 
served in the COSMOS (jllbert et all 120091 ) field; this cat- 
alogue goes to an i-band depth of 25 mag, but we limit 
our detection depth to an i magnitude of 22 to allow for 
magnifications of ~ 15 by strong lensing. Both the optical 
and millimetre depths are conservative in light of upcoming 
telescopes. However, these represent the limit to which we 
trust the models; beyond these depths a currently unknown 
population of dim objects may exist, although if the shape 
of the redshift distributions do not vary significantly with 
flux limit then the cosmological conclusions are indepen- 
dent of survey depth. Given a population of known lenses 
we evaluate the probability of being a double source plane 
lens, P(D\L): for a [z\, cry] pair we calculate the lensing 
cross section and magnificati on tensor as functions of red- 
shift, and apply these to the lBethermin et al.l (|201ll ) model 
or COSMOS catalogue. This me t hod i s broadly similar to 
the method in lOguri fc Marshall! (|2010n . For axi-symmetric 
lenses our method is equivalent to solving for each lens the 
integral 



P(D\L) = 



roo re 

/ dz 

J z, JO 



dSa CIOBB (z, zi,a v ) 



d 2 N 
dzdS 



W(S,Su. 



(10) 
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where ^ z ^ s is the differential number distribution of second 
sources per unit solid angle with respect to redshift and 
flux, across is the lensing angular cross section, and W is a 
function of the flux (S) and the survey flux limit (Sum); W 
encodes the magnification bias. In the case of an SIS, 



3 (z,Zi,av) = it \t 



,SISf 
B 0, 



z \"v) 



(11) 



and 



[ 1 if S > Slim/2 

«'<«»"•= {(».)' ifs<Sll „ /2 < 12 > 

with Equation 1101 gives the probability for finding a second 
source behind each lens. The probability distribution func- 
tion (PDF) for finding a source at redshift 2 S 2 is the deriva- 
tive of Equation \W\ with respect to z. We then generate 
the second source redshift by drawing uniformly from the 
PDFs to generate a population of [z\, z s i, z S 2] systems, with 
a weight given by Pi(D\L). Given the weighted population, 
we sample 1000 sets of n systems from the population, and 
forecast the limits on w for each set. These 1000 limits are 
combined to produce Figures [&I 4101 the histograms show the 
distribution of upper and lower 68% confidence limits. The 
medians are taken as the expected limits and the 16th and 
84th percentile are taken as the errors upon the expected 

limits (68% CL}. 

Using the lBethermin et al.l (|201lh model and the veloc- 
ity dispersion function of a modified Sch echter function fit 
to the Sloan Digital Sky Survey (SP SS; lYork et al.ll200oT ) 
data bv lChoi. Park, fc Vogelevl (J2007I 1. 



dn = 4>* I — 



cxp 



g da 



(13) 



where (/>» = 8.0 x 10~ 3 /i 3 Mpc~ 3 , a* = 161 kms -1 , 
a — 2.32 and f3 — 2.67 for galaxy scale lenseqj, we calcu- 
late that typically 1 in every 250 red elliptical s will be a lens 
for sou rces with S>0.3 mjy at 1.1 mm. Whilst IGavazzi et al.l 
1)20081 ) showed that lenses represent a population with bias 
towards higher mass, one can take their calculations of this 
bias to estimate that roughly 1% of typical galaxy lenses will 
yield a second source with S>0.3 mjy at 1.1 mm. As such we 
find that a blind search of galaxy scale lenses is unlikely to 
yield double source plane systems unless deeper sensitivities 
are probed. However, by looking at surveys that preferen- 
tially select more massive systems (groups, or the rare very 
massive elliticals) we ca n increase that rate. A s an example, 
the Cosmic Horseshoe (jBelokurov et al.l 120071 ) has a ~30% 
chance of a S>0.3 mjy (at 1.1 mm) object lying behind it, 
whilst a [z — 0.2, a = 350kms _1 ] elliptical will lens a 
S>0.3 mjy source ~6% of the time. 



4.1 Follow up of a photometrically selected wide 
arc catalogue 

The Cambridge and Sloan survey of wide arcs in the sky 
(CASSOWARY) targets wide separation blue arcs around, 
or multiple blue companions to, luminous red galaxies 



wc do not allow this function to evolve with redshift 



Model CASSOWARY 
Actual CASSOWARY 




0.4 0.6 

lens redshift 

Figure 7. The lens population of our CASSOWARY-like cata- 
logue (shown in grey), compared to the known redshifts of 106 
CASSOWARY lenses (shown hatched; not all of these are avail- 
able in the currently published literature). The simple model de- 
tailed in Section |4TI broadly replicates the known CASSOWARY 
lens redshift distribution. 



(LRGs) in the SDSS imaging catalogue. This method pref- 
erentially selects the most massive lenses, but it also often 
selects complicated group-scale lenses that can be difficult 
to model; this is the very problem that makes the double 
sourc e plane method difficult for clusters (but see I Julio et al.l 
[201fJ). For small groups it may be more tractable and the 
large, bright images that are typical in CASSOWARY pro- 
vide a lot of positional information for lens models to fit 
and the larger lens masses that generate these wide arcs 
greatly increase the prospects of finding a second source. 
The wide arcs tend to be at relatively high redshifts, up 
to about z = 2, although there are also some arcs with 
z < 1. For the high z s \ sources, a useful follow up strat- 
egy is to look for objects intermediate between known lens 
and source (probably in the optical), whilst the population 
of low redshift wide arcs should be followed up to find a 
high redshift background source (probably in the millime- 
tre). Selecting the optimal transition point depends on the 
position of the lens and known source, but for simplicity we 
make forecasts for the constraints possible from observing 
the whole catalogue in the optical and in the millimetre. 

Much of the CASSOWARY deep optical imaging fol- 
lowup is already being conducted to study the lens and 
source (i.e., for non-cosmographical reasons) and this is 
likely to continue. Several of the lenses are known to contain 
two or more unique optical sources, though spectra have not 
always been obtained for the sources. The CASSOWARY 
catalogue is still undergoing spectroscopic follow up, so there 
are currently few systems with known <j\ and z s i; this pre- 
cludes us from estimating the CASSOWARY selection func- 
tion simply by taking the known systems [zi,Zsi]. Instead, 
we generate a mock CASSOWARY-like system of group- 
scale lenses. We assume a constant comoving density of po- 
tential lenses and take the luminosity distribution for the 
central galaxy as t he red ellipticals luminosity function of 
iDriver et al] J2007r ) 



dn 

dM 



oc ^■MM-M, ) y+' exp j_ (y.4(M,-M)\j (14) 

where M is the absolute B-band magnitude, M, = — 19.02 + 
51og 10 (/i) and a = —0.26. The lensing cross section of each 
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-1.4 -1.2 -1.0 

Upper and lower la limits on w 

Figure 8. The upper and lower bounds on the 68 percent con- 
fidence region, obtained using 1000 selections of 10 (red), 20 
(hatched white) and 40 (green) double source systems. The dou- 
ble source systems were generated using the CASSOWARY-like 
catalogue and the COSMOS photometric redshift catalogue down 
to «i; m = 22. There are no green bars concealed behind the red 
bars. 
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Figure 9. As Figure \E\ but with 20 double source plane sys- 
tems and different uncertainties on rj: a(r])/r] = 1% in green, 
t7 ( r l)/ r l = 2% in hatched white, a(r])/r] = 4% in red, 
cr ( r ?)/ r ? = 10% i n blue. The x scale is the same as in Figure [8] 



potential lens is then projected onto the COSMOS popula- 
tion in the same way as Equation [TO] except we are evaluat- 
ing the probability of being a lens P%{L) in place of Pi(D\L). 
The probability of being a CASSOWARY-like lens requires 
us to make selection cuts on our model lens-source sam- 
ple and a system is only included if it passes the following 
cuts: lens apparent magnitude, r < 21.5; Einstein diameter 
greater than 3"; a limiting total arc magnitude of r < 22; arc 
colour of g — r < 0.5; lens absolute magnitude Mb < -21.5. 
The last cut is made because typical LRG selections exclude 
low lumin osity objects, as is seen clearly in the luminosity 
function of lEisenstein et al.l i|200ll ) . Once this procedure has 
been applied, we are able to create a weighted population 
of <j\ and z s i for each z\. Our catalogue broadly replicates 
the known lens redshift of 106 CASSOWARY lenses (and 
lens candidates), as shown in Figure [7] Since the purpose of 
generating this population is to estimate future cosmological 
constraints rather than to perfectly replicate any particular 
current survey, we are unconcerned by the simplicity of this 



Table 1. Expected la upper and lower bounds on to, given n 
CASSOWARY double lens systems. The lens and first sources 
are generated using the CASSOWARY-like model and the sec- 
ond sourcescomefromthe deep optical catalogue (i < 22), and 
the lBethermin et al.l J201 J) millimetre model (S > 0.3 mjy). The 
Section marked "Combination" has second sources taken from 
both the second source populations, with a weighting of approxi- 
mately three optical second sources per millimetre second source. 
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CASSOWARY-like model. With this mock catalogue of [z\, 
a"i] we can generate second sources using the method out- 
lined above, the weights for each double source plane system 
are now given by Pi(L) x Pi(D\L) if the lens and first source 
pass all the selection cuts, and is zero otherwise. 

If the lenses are assumed to have an effective velocity 
dispersion of (a 4 ) 3 Ri 350 kms -1 , we find 7% to also be lens- 
ing an S>0.3 mjy millimetre source. For the CASSOWARY 
objects already observed, 3 of 15 have so far been discovered 
to have a second arc. We assume therefore that 7 millime- 
tre sources and 20 optical sources would be yielded by a 
fullscale follow-up campaign. We make forecasts for the cos- 
mographic constraints possible with 5, 7, and 15 millimetre 
sources and with 10, 20, and 40 optical double source sys- 
tems. Figure [8] shows the likelihood distribution of the up- 
per and lower limit on w with 10, 20, and 40 optical second 
sources (a(ri)/ri — 1%). Finally we combine the datasets 
(still assuming a(rj)/r] = 1%) to draw 15, 30, and 50 dou- 
ble source systems. Because of the complexity of group-scale 
structures it may prove impossible to measure r\ to 1% un- 
certainty; Figure [9] shows the likelihood distributions for 20 
double optical sources if a 1%, 2%, 4% or 10% uncertainty 
on r\ can be achieved; the results are similar to Figure [3] 
in that the constraints on w depend significantly upon o{rj), 
but given 20 systems even 4% uncertainties on 77 can provide 
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Table 2. Expected la uncertainty on w, given n SLACS dou- 
ble lens systems, detected at a flux limit of S. In all cases 
<r(r,)/r, = 1%. 

S = 1 mjy 



„ _n 72 n+ a030 
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S = 0.3 mjy 
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meaningful constraints on w. If cr(r])/rj is as high as 10% a 
population of many more than 20 double source plane sys- 
tems would be needed to make an important contribution to 
cosmography. All of the results for this CASSOWARY-like 
model are shown in Table [TJ 
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Figure 10. The upper and lower bounds on the 68 percent confi- 
dence region, obtained using 1000 selections of six double source 
systems. The double source systems were generated using the 
SLACS catalogue and the Bethermin et al. model for 1.1 mm 
continuum observations the grey bars are for sources selected at 
S>0.3 mjy and the red bars are from sources selected at S>1 mjy. 
The constraining power degrades by ~1 percentage point, for 
S>0.3 mjy compared to the S>1 mjy results. 



4.2 Follow up of a SLACS type catalogue 

The Sloan Lens ACS (SLACS) Survey selected strong lens- 
ing candidates using the SDSS spectroscopic survey. Specif- 
ically targeting bright early type galaxies, SLACS takes 
advantage of the 3" diameter spectroscopic fibres used by 
the SDSS and looks for emission line features that sug- 
;est a higher redsh i ft component than the target galaxy 
Bolton et all I2004J . 120061 . I2008J : lAuger et al] I2009D ; these 
targets are then followed up using ACS and WFCP2 on the 
Hubble Space Telescop e (HST). The sele ction function of 
SLACS is complicated i|Dobler et al.ll2008T ). but large Azi_ s i 
systems are selected against since the emission lines are of- 
ten redshifted out of the bandpass, and the finite diameter of 
the fibre requires small Azi_ s i for the most massive lenses. 
As Section [2] showed, bias towards low Azi_ s i is desirable 
for finding optimal double source systems. The selection of 
lower mass objects (compared to group-scale lenses, for ex- 
ample) decreases the prospect of finding a second source 
behind the SLACS system because the lensing cross section 
scales with mass squared; however the ideal follow up can- 
didates (small Azi_ s , high mass) still show up in the SLACS 
sample-because of th eir small Einstein radius-despite their 
high mass. Using the IBethermin et al.l (|20111 ) model, a full 
millimetre follow up of the current SLACS catalogue down 
to 1 mjy should produce ~1.5 new double source systems, 
whilst observing to 0.3 mjy would yield ~3.0. Since the cat- 
alogue is growing, and this is a conservative flux limit given 
the observing capability of ALMA, we forecast the expected 
constraints on w given by 1, 2, and 6 double source systems 
found with S>1 mjy and 3, 6, and 12 systems found with 
S>0.3 mjy. The results are summarised in Table[2] FigurefTOl 
shows the likelihood distribution of the upper and lower lim- 
its on w given six double source systems with S>0.3 mjy and 
S>1 mjy. Figure [TT1 shows the w-Qm plane for a typical se- 
lection of 6 SLACS-type S>0.3 mjy double sources; the six 
"typical" systems are listed in Table [3] 

Directly comparing the two instances of six double 



Table 3. The typical set of 6 SLACS-type lenses. These arc the 
6 systems that gave the median constraints on w for WMAP+6 
double source plane lenses in Table [2] One SLACS system has 
reappeared (z\ = 0.227, z B = 0.931), but such repeats are not 
excluded by the weighted selection of lens-source systems; neither 
Systems 2 or 4 are likely to play a significant part in constraining 
w due to their large Azi_ sl . For cosmography, System 5 is the 
best system by some margin, although systems 1 and 3 are also 
significant when constraining evolving models. 



System 


1 


2 


3 


4 


5 


6 


zi 


0.440 


0.227 


0.195 


0.227 


0.194 


0.111 


Zsl 


1.192 


0.783 


0.632 


0.931 


0.446 


0.316 


Zs2 


3.859 


0.931 


2.724 


1.766 


2.058 


0.463 



source lenses shows some degrading of the constraint on w as 
the flux limit is lowered. This is due to th e large number of 
low flu x, low redshift objects predicted by IBethermin et al] 
i|2011T l which are individually less useful for cosmography 
than the higher-redshift systems. A shallow survey of lots 
of lenses would select fewer of these low redshift millimetre 
sources (and thus be optimal for cosmography), but since 
there are very few known (simple) high mass lenses, such a 
survey is currently impossible. Despite the fact that SLACS 
provides much better cosmography per double source plane 
system than the CASSOWARY-type, our results are dom- 
inated by the few SLACS objects with very high velocity 
dispersion and close first source: the lower velocity disper- 
sion objects are unlikely to be double source plane systems 
and there are no high velocity dispersion systems with in- 
termediate or large lens-source separation. Indeed, many 
high mass lenses will be missed by SLACS due to the source 
arcs falling outside the spectroscopic fibre. These high mass 
systems will be crucial if a large number of double source 
systems are to be discovered efficiently, so follow up obser- 
vations of all the (simply-modeled) high mass systems in 
strong lens catalogues should be conducted. 
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Figure 11. The w-(!m plane for a typical selection of 6 SLACS- 
type double source plane systems (red), with WMAP (grey) and 
the combination of both (black). The six systems are listed in Ta- 
ble [3] The confidence regions have shrunk significantly compared 
to one system (Figure [2]l and are still orthogonal to the WMAP 
constraints. 



5 BEYOND WCVM 

Up to this point we ha ve only investigated models with con - 
stant equation of state. ICaldwell. Dave, fe Steinhardtl (|1998l ) 
first suggested the possibility of an evolving equation of state 
in the context of a quintessence model. We now investigate 
the power of double source plane lenses to constrain w(z) 
models. We adopt the parametrisation 



w(z) — wo+w a (l — a) 



— w + w a 



1 + 2 



(15) 



as int roduced in IChevallier fc Polarskil i|200lr ) and iLinderl 
(2003). Our fiducial cosmology is a special case of this 
parametrisation with rao = — 1 and w a = 0. This 
parametrisation allows us to solve the integral in Equation 
[()] analytically: 



JM 



= (1 + *) 



3(l-\--w +w a ) 



exp 



-3w a z 
l + z 



(16) 



Taking the 6 systems in Table [3] cr{rf)/r) = 1%, and uni- 
form priors [0 < Qm ^ 1; — 3 ^ Wa ^ — i; — 3 < w a ^ +3], 
our MCMC forecast finds the confidence intervals shown in 
Figure 1121 In high £7m universes, the dark energy equation 
of state is poorly constrained by measuring r\. Marginal- 
ising over all Qm leaves very weak constraints in the wq— 
w a plane because of the high Q,m samples. With a prior 
on Qm < 0.5 the wo-w a plane shows a boomerang shaped 
degeneracy, which arises from taking the ratio of angular di- 
ameter distances. A universe with Wo + W a > has no era of 
matter domination and is strongly excluded by CMB con- 
straints; coincidentally this is the region for double source 
plane constraints that is allowed by taking angular diameter 



s o 




Figure 12. Forecast parameter constraints for the w(z) = too + 
w a (l — a) model. Constraints are from the six systems listed in 
Tablc|3]assuming <j(r])/ri = 1%, and uniform priors [0 Sj Hm ^ 1; 
—3 ^ too ^ 0; —3 ^ w a ^ +3]. The plot shows 1, 2 and 3 a 
regions in the too — w a plane, marginalised over SIm, the samples 
are coloured by fijyi ; a stronger prior on S!m is needed to improve 
the constraint on too and w a . 



distance ratios. The six double source plane systems alone 



give wo = 



1.177 



+0.460 



= 0.485l2;353. Combined 



with a CMB prior - forecast for Planck - the 1 a constraints 



are wq 



-0.802 



+0.495 „„ _ _ n rkk+1.228 

-1.175- 

For the Planck MCMC forecasts, we use the 
seven frequency channels that have polarisation sensi- 
tivity^ with noise levels given in the Planck Bluebook 



l|The Planck Collaborationll2006f ). In addition we include a 
reconstruction of the CMB weak lensing potential, with sta- 
tistical n oise calculated using t he optimal quadratic esti- 
mator of lOkamoto fc Hul (J2003J). We use information from 
the spectra up to a maximum angular multipole of 2500, 
and assume that 65% of the sky is observed. We sam- 
pled the likelihood using the publicly-available CosmoMC 
package (ILewis fc B ridle 2002]) , with fiducial spectra calcu- 
lated with CAMB (|Lewis et alj|2000h . Evolving dark energy 
models wer e impl emented with the PPF CAMB module of 
iFang et all (120081 ). In Figure [13] we show the Planck and 
6 double source plane system constraints for each parame- 
ter pair in the flat w a model, marginalized over the third 
parameter. 

It is clear from Figures \T%\ and [TH] that wo and w a are 
degenerate; as such we use the figure o f meri t (FOM) as 
defined in Equation 6 of lMortonson et al.l (J2010|) to quantify 
the constraining power: 



FOM = 



6.177T 

^95 ' 



(17) 



where Ag 5 is the area enclosed by the 95% confidence in- 
terval in the wo~w a plane marginalised over all other pa- 
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Figure 13. Forecast cornerplot for the w(z) = wo + u> a (l — a) model. On diagonal subplots shows the marginalised likelihood for each 
parameter. Off-diagonal subplots show the two-dimensional 1 and 2 a confidence contours, marginalised over the third parameter. The 
grey contours are forecast for Planck (CMB lensing and ISW effects have been included). Olive shows the constraints from 6 double 
source plane systems in Table \3\ assuming <j(rj)/rj = 1%. Black contours are forecast for the combination of Planck with the 6 double 
source plane systems. 



rameters. Using this definition the FOM for Planck plus the 
6 lenses in Table [3] is 14.2. Using WMAP+6 systems low- 
ers t he figure of merit to 7.6 . When asserting spatial flat- 
ness, iMortonson et all (J2010I ) find FOM = 15 for WMAP 
plus the Union supernovae compilation (jKowalski et al.l 
2008), combined w i th th e baryon acoustic oscillation data 
o nEisenstein et al.l l|2005l ) and an Ho prior. 

Figure [T4l shows the FOM contour plot for Planck plus 
one double source plane system with a(rj)/r] — 1%; we 



find that for wq , w a models the figure of merit increases 
as z S 2 increases and as A2i_ a i decreases. We also find that 
the wo — w a degeneracy increases (away from vertical in 
Figure I13p with increasing lens redshift, conversely below 
2i ~ 0.2 the degeneracy is almost broken, but the FOM 
rapidly degrades as z\ decreases. Since the degeneracy be- 
tween parameters changes with the redshift of lens and 
sources, the FOM is likely to grow quickly with sample size. 
To strongly constrain evolving models it will be important to 
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Figure 14. Figure of merit contour plots for the w(z) = wo + «i a (l — a) model, with FOM = 6.I771-/A95. Constraints are for a single 
double source plane lens assuming cr(r])/r] = 1%, combined with a forecast Planck dataset. Left Lens redshift z\ = 0.35. Right Lens 
redshift z\ = 0.55. A large second source redshift and small separation between lens and first source gives the best FOM. 



have a good coverage of lens re dshifts. The ongoing BELLS 
survey ( Brownstein et al.l 12012 ) uses very similar selection 
methods to SLACS, but targets higher redshift potential 
lenses; their 25 definite and 11 probable galaxy-galaxy sys- 
tems have lens redshifts 0.4 < z < 0.7. These higher redshift 
systems should be especially useful for incressing the double 
source plane FOM - if a second source can be found behind 
any of them. 

Throughout this work we have asserted a spatially flat 
universe, but allowing f2j. to be free weakens the constraints 
on w. Without asserting spatial flatness we find that WMAP 
plus the 6 typical systems give the following results for a 
constant equation of state: w = — Filing. For non- 
flat evolving equation of state models our FOM for Planck 
plus 6 lenses degrades slightly to 12.3. Double source plane 
systems place no significant additional constraints on Slfe 
once the Planck results are included due to the excellent 
precision with which CMB lensing constrains flatness. 



6 DISCUSSION AND CONCLUSION 

Strongly lensed double source systems and the CMB are not 
the only datasets to constrain w; Supernovae (SNe), Baryon 
Acoustic Oscillations (BAO), galaxy clusters, and time de- 
lays in strongly lensed time variable sources (At) also pro- 
vide constraints on w. Whilst only the strongly lensed dou- 
ble source method is independent of the Hubble constant, 
WMAP+BAO+Ai can still give strong constraints on w 
when combined with a prior on h obtained using local su- 
pernovae (we have not incl uded high redshift sup ernovae or 
galaxy cluster constraints). iKomatsu et al.l 1)20111 ) used this 
combination of methods, finding w = — 1.08 ± 0.13. In the 
future the constraints from strongly lensed double source 
plane systems can easily be combined with other datasets 
to give an improved constraint on w. The main advantages 
of this are the breaking of degeneracies and, perhaps more 
importantly, the investigation of systematics inherently af- 
forded by using different measurements. Figure 1151 shows 
the w-flM plane for a typical set of 6 SLACS type double 
sources (listed in Table [3| given a uniform flat prior on w 
and Om and a forecast combination of WMAP+BAO+Ai+6 
double source plane lenses. In this case we recover a con- 
straint of w — — 1.041^9 , for our fiducial cosmology of 
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6 Double source plane lenses 
WMAP+BAO+At+6 Doubles 




WMAP+BAO+A£+6 Double source plane systems 
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Figure 15. Top The ui-Om plane for a typical selection of 6 
SLACS-type double source plane syst ems (red), with WMAP 
(grey) and the BAO constraints of Eisens tein et al.l (2005) 
(dashed blue). The combination of W MAP+BAO+Af. and the 
6 SLACS lenses is shown in black. The lSuvu et al.l J2010I ) results 
are not shown in this figure, but are included in the black contour. 
Bottom The marginalised distribution of w for WMAP+BAO+ At 
(hatched) and the improvement gained when adding 6 SLACS- 
type double source plane systems (grey). 
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w — — 1. Setting the fiducial cosmology to w = — 1.08 
yields w = — 1.08JIn'nQ- These represent a 30% improvement 
on the iJKomatsu et al.l l201lh result, although the At and 
BAO constraints will improve significantly as their datasets 
grows and if h is better constrained; this method might be 
slightly less competitive by the time 6 double source plane 
systems have been found, but a significant strength of mul- 
tiple source plane lenses is as a constraint on w independent 
of any measurement of h. 

Most w ,w a forecasts have used a Fisher for- 
malism to estimate f u ture parameter con straints (e.g 
ICoe fc Moustakas! 120091 ; IShapiro et all l2010h . Whilst the 
Fisher forecasts are not directly comparable to our MCMC 
forecasts it seems that the wo,w a degeneracy for double 
source plane lenses will be in a similar direction to BAO, 
SNe and weak lensing. For all these probes low wo corre- 
lates with positive w a - However, for very low lens redshifts 
(z\ w 0.2), the double source plane constraints on wo become 
almost independent of w a , so there is scope for complemen- 
tarity, but at the expense of the constraining power of each 
lens. 

The problem of finding second sources will not be triv- 
ial, and the probability of finding a second source scales with 
mass squared. As such a dedicated search strategy should 
look to find second sources behind high mass lenses. Cur- 
rently there are a handful of good candidate systems, and 
a few double source systems should be uncovered with a 
moderate investment of telescope time. With a larger popu- 
lation of known lenses, the number of very good candidates 
for such a search will increase, making searches for com- 
pound lenses more efficient and more useful; not all double 
source systems will be particularly useful for cosmography, 
mostly due to having both sources far beyond the lens, but 
a selective search strategy can minimise the time spent on 
these objects. One wants to draw the second source from a 
different population to the first to make z s i ~ z S 2 less likely; 
this can either be achieved by looking in a different wave- 
length or to significantly greater depth. It may also not be 
possible to measure n to 1% for all objects, especially af- 
ter accounting for lensing by the second source. However, 
the CASSOWARY results for 20 lenses with a 4% uncer- 
tainty on r\ show that given a large sample, w can still be 
constrained with double source plane lens systems. Further 
work will evaluate the systematics of measuring 77, paying 
particular attention to the mass model of the lens and lens- 
ing by the intermediate source or other perturbers along the 
line of sight. 

(i) Can double source plane lenses be used for cosmogra- 
phy? 

Double source plane lenses can be used for cosmography. 
We have shown that - uniquely amongst cosmological 
probes - it is possible to constrain the equation of state 
with double source plane lenses, independent of the Hubble 
parameter. 

(ii) What is an optimal configuration of lens and source 
redshifts for cosmography? 

By resampling the WMAP team's MCMC chains we have 
shown that, for cosmography, the key features are a low 
separation lens-source pair and a high redshift background 
source. We found cosmography to be broadly insensitive to 



lens redshift, so long as it is greater than about 0.2. 

(iii) How well could cosmography be constrained with a 
realistic population of double source plane lenses? 

Given our knowledge of what represents a good system for 
cosmology, we postulate a search strategy to look for double 
source plane systems; ultimately such a follow up strategy 
will depend on the single source plane systems (which will 
li kely increase in number ), but using the population model 
of lBethermin et al.l (|2011f ) and the photometric redshift cat- 
alogue of COSMOS, our model allowed us to make forecasts 
for a population of double source plane systems discovered 
in a variety of different ways. From this population we have 
learnt that - at a{r\)jr) = 1% - only a handful of dou- 
ble source plane systems are needed to be competitive with 
current constraints on w and that 6 systems can provide 
meaningful constraints for wo,w a models. The JIm-w de- 
generacy of this method is almost orthogonal to CMB and 
BAO datasets, making the double source plane method com- 
plimentary to existing work. 
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